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We investigate the effect of gravitationai iensing, produced by linear density per- 
turbations, for anisotropies of the Cosmic Background Radiation (CBR) on scales 
of arcminutes. In calculations, a flat universe (0 = 1) and the Harrison- Zel'dovich 
spectrum (n = 1) are assumed. The numerical results show that on scales of a few 
arcminutes, gravitational Iensing produces only negligible anisotropies in the tem- 
perature of the CBR. Our conclusion disagrees with that of Cayon et al. who argue 
that the amplification of AT/T on scales < 3' may even be larger than 100%. 
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I. INTRODUCTION 

Many authors have computed the effect of gravitational Iensing for anisotropies of the 
CBR (Blanchard & Schneider 1987; Kashlinsky 1988; Cole & Efstathiou 1989; Tomita & 
Watanabe 1989; Linder 1990; Watanabe & Tomita 1991; Feng & Liu. 1992; Cayon et al. 
1993). Unfortunately, however, their conclusions are controversial. Roughly speaking, there 
exist three different kinds of conclusion so far. The first is that gravitational Iensing effects 
strongly erase fluctuations of the CBR on scales of a few arcminutes (Kashlinsky 1988); the 
second is that an appreciable, even strong, amplification of AT/T is possible (Sasaki 1989; 
Linder 1990; Cayon et al. 1993); the last is that gravitational lens effects on the CBR are 
negligible (Cole & Efstathiou 1989; Tomita & Watanabe 1989). 

Recently, Cayon et al. presented calculations of the gravitational Iensing effects, pro- 
duced by linear density fluctuations on the CBR and got an interesting result. Their work 
implied that there should be an appreciable amplification, of the order of 20%, for AT/T in 
present experiments on the scales of several arcminutes, whereas previous work including the 
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effect of nonlinear density fluctuations found a negligible amplification (Cole & Efstathiou 
1989). This implies that the gravitational lensing effect on anisotropics of the CBR due 
to linear density fluctuations overwhelms that due to nonlinear clustering. This result is 
surprising and difficult to interpret. 

In this paper, we use a new formalism to calculate the effect of gravitational lensing 
on anisotropics of the CBR produced by linear density fluctuations. Our calculations show 
that the gravitational lensing effect on scales of a few arcminutes is essentially negligible in 
contrast with the results of Cayon et al. 

Within the geometrical optics approximation, a gravitational field is equivalent to an 
optical medium with a refractive index different from unity, and the deflection of light may 
be interpreted in terms of the refractive index and its spatial variation. The amplitude and 
phase fluctuations, produced by a random gravitational field, may then be calculated using 
the usual methods of random medium optics (Fang 1982). 

In Sec. II, we establish the basic equations for a wave scattered by a gravitational potential 
and for the amplitude of the scattered wave. In Sec. Ill, we present the explicit formulae 
for calculating the anisotropics of the CBR. Numerical results and brief conclusions are 
summarized in Sec. IV. 



Under some reasonable assumptions (Cole & Efstathiou 1989), one can prove that there 
exists a gauge such that the metric perturbations are characterized by a single potential 



where a(t) is the scale factor of the universe. The relationship between the gravitational 
potential <f> and the matter density perturbation Sp is with 8ttG = c = 1 



II. THE BASIC METHOD 



<t>{t,-x) < 1 





A(f> = -a 2 
* 2 
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In the background universe for the linear perturbation (which we assume to have fi = 1), 
<f> is time-independent. 

For the sake of simplicity in writing the Maxwell equations in a perturbed metric, we 



The speed of light in the metric of Eq.(3) is set to unity. A fundamental problem for treating 
light propagating in an inhomogeneous universe is that of how to describe the gravitational 
lensing effect. According to geometrical optics, a gravitational field is equivalent to a medium 
with a refractive index different from unity and the deflection of light may be interpreted in 
terms of the refractive index and its spatial variation. As a simple version of this analogue, 
consider the spatial part of the photon four-momentum = to be directed along the 
x 3 axis; one then has the following geodesic equation: 

dV = -2^dx 3 , i = l,2. (4) 

dx % 

On the other hand, the change in the direction of a light ray propagating in an inhomoge- 
neous medium with refractive index n = 1 + n 1 is 

dk' = —dx 3 , 1 = 1,2. (5) 

dx l J 

By comparing of Eqs.(4) and (5), it is obvious that the effect of the perturbed gravitational 
field is equivalent to a change in the refractive index, i.e. 



Therefore, considering the equivalence of the two descriptions of light propagating in an 
inhomogeneous matter distribution and in a medium with inhomogeneous refractive index, 
the Maxwell equations may be written in the following form 



use the conformal time r 



3v / 3^ 1/ ' 3 instead of t. Then 





ri = -2<p 



(6) 




(7) 
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^ = -(l + 2^)VxE. (8) 

Formally, adopting Hanni's approach (Hanni 1977), the Maxwell equations in an universe 
with an inhomogeneous matter distribution may be obtained in more rigorous way. 

The propagation of an electromagnetic wave in a random inhomogeneous medium is 
accompanied by a number of fluctuation phenomena including polarization, fluctuation of 
phase and fluctuation of amplitude. In the present paper, we concentrate on fluctuation of 
amplitude. 

In a flat universe, because the gravitational potential is time-independent, we can derive 
a solution of the Maxwell equations (7) and (8) representing monochromatic waves with 
fixed frequency u. In other words, we can assume that the electric and magnetic fields have 
the form i?e(Ee~ 8lx ' T ) and i?e(He" !WT ) respectively. In the case of A = 2ttlo~ 1 <C / , (7o is the 
typical length-scale of inhomogeneity for matter in the universe), keeping only first order 
terms, the wave equation can be simplified as (Feng and Liu 1992) 

(V 2 +oj 2 )E = -4w 2 (f>E. (9) 

To solve the equation of wave scattering, in the case of weak fluctuations, we apply the 
Born approximation to expand E in the series 

E = E + E 1 + E 2 + • • • (10) 

with Eq 3> Ei ^> E 2 3> • • •• Putting Eq.(10) into Eq.(9), we get two first-order equations: 

(v 2 + ^ 2 )£o = o, (li) 

( V 2 +^ 2 )£i = -4cu 2 ^o- (12) 

If the CBR is perfectly uniformly distributed, there is no net gravitational lensing effect 
on anisotropics of the CBR. This is a well-known result in geometrical optics and will be 
confirmed below from the point of view of wave scattering. In order to study the gravitational 



lensing effect on anisotropics of the CBR, we focus on how the fluctuation part of the CBR 
is affected. To do so, we separate E into a homogeneous part E h and a fluctuation part 

E = E oh + E f (13) 

with 

Eoh =< E oh >; E f = E — < E f > . (14) 

Here, the averaging is done over all observation directions. As Eoh and E j propagate 
in an inhomogeneous universe, they interact with the gravitational potential and produce 
scattered waves E\h and E\j respectively. Eoh and E j and their scattering terms satisfy 
following equations 

(V 2 +w 2 )£ofc = 0, (15) 

(xj 2 + to 2 )E lh = -4w 2 <f>E oh , (16) 

( V 2 +^ 2 )^o/ = 0, (17) 

(V 2 +u 2 )E lf = -4u 2 ^E of . (18) 

After scattering, the outgoing wave is 

E = E + E 1 = E h + E of + E lh + E lf . (19) 

Because the last scattering surface is far away from us, without loss of generality, we let the 
incident waves have the plane wave forms, Eoh = Aoh^'* and E j = A /e 8lx " x . The scattered 
waves can then be easily expressed, using Green's function method, in the following forms 

= — / 44x)<£(x) ^ 3 x, (20) 

7T J X 



£i/(xi) = - / 

7T J 



A 0/ (x)^(x) 



,i(iOX — LU-X.) 

J3 



<fx. 



(21) 



x 



When Eqs.(20) and (21) are evaluated, it is sufficient to include only the contribution due 
to waves scattered through angles not exceeding 9 = <C 1. In other words, the integration 
can be confined to the part of space which lies within the cone C(dVl), 9 < j^, where 9 is 
the angle between the direction of observation and the direction of the scattering element, 
and l is the typical inhomogeneity scale in the universe. In fact, the integration functions 
of Eqs.(20) and (21) oscillate rapidly outside the cone C(o?fi), so that for a sufficiently 
smooth variation of </>(x), integration over the region external to cone C(dft) only provides 
a negligible contribution. This cone C(dft) is much smaller than the angular scale of CBR 
inhomogeneities and so we have 



A)/i(xi) is the amplitude of the uniform part of the CBR and A /(xi) is that of the fluc- 
tuation part; as a direct consequence, Eih{x-i) is independent of observation direction, but 
£i/(xi) does depend on the observation direction. 



In the Rayleigh- Jeans part of blockbody radiation spectrum, T oc I (where T and I are 
temperature and intensity of the CBR respectively) and I oc A 2 (where A is the amplitude 
of the electric vector E). We will not distinguish A and E in the following except when 
necessary. We then have 






(23) 



III. ANGULAR CORRELATION FUNCTION 




(24) 



Putting Eq.(19) into Eq.(24) and using the definitions of Eoh and Eqj, we have 
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It is obvious from Eqs.(22) and (23) that E\h ~ (Eih) = and (E\f) = 0. Without loss of 
generality, we let Eoh = 1 and then the above equation simplifies to 

S -^ = E of + E lf . (26) 

Considering the above relations, the angular correlation function of anisotropy of tempera- 
ture of the CBR may be obtained 

C(a) = { mMJ_ipiy (27) 

where cos(a) = Xi -x 2 and the averaging is done over all observation directions. Substituting 
Eq.(26) into Eq.(27), it follows that 

C(a) = (| 2£ 0/ (x 1 ) | 2 ) + 2(| 4£ 0/ (xi)£i/(xi) I) + (I 2E 1/ (x 1 ) | 2 ). (28) 

In the right-hand side of Eq.(28), the first term is the angular correlation function of the 
primordial fluctuation background; the second describes the interaction between the gravita- 
tional lensing effect and the perturbed part of the primordial background (this is called the 
angular cross correlation function C c (a), and determines the lensing effects on the anisotropy 
of the primordial CBR); the last is a higher order term and may reasonably be omitted. The 
major purpose of this paper is to determine the angular cross correlation function C c (a). 
We Fourier decompose the gravitational potential, obtaining 

<K X ) = (3^3 / A k e 8k - X ^ 3 k. (29) 

Based on linear perturbation theory and assuming a Gaussian random fluctuation field, we 
have 

A k = ^H 2 k~ 2 S k , (30) 
(hM=\6 kl f £ 3 (k!-k 2 ), (31) 
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(A kl A k2 ) = ^k~ 4 | 6 kl | 2 5(k! - k 2 ), (32) 

(«5 kl A k2 ) = hilk- 2 | 6 kl | 2 5(k! - k 2 ), (33) 

where | ^ | 2 and | | 2 are the spectra of perturbations of the density and gravitational 
potential respectively. 

Roughly speaking, for adiabatic perturbations, the temperature fluctuation of the CBR 
on scales of a few arcminutes is 6T/T = (l/3)6p/p before the recombination era (Silk 
f 967). What are the amplitude and shape of the temperature fluctuations after decoupling? 
Of course, this depends very much on the assumed recombination history. For the sake 
of simplicity, we assume that the time-scale of recombination is extremely short so that 
ST/T = (l/3)6p/ p remains with the previous amplitude and shape. This assumption is rea- 
sonable for our purposes since we are dealing with wave propagation and its interaction with 
gravitational lensing between the last scattering surface and the observer. Our assumption, 
in fact, just means choosing a convenient initial condition. 

Fourier decomposing the fluctuation part of the CBR and its scattering term, we obtain 

£o/(xi) = ^-y 2 Js^fk, (34) 

2£ 1J (x 1 ) = | J I^Ak.e^^e^^cos^-cu-x)^ 3 ^^, (35) 

where the vector Xi points to the last scattering surface and has length 2i7 ( 7 1 ; the unit 
vector Xi represents the direction of observation. 

For mathematical convenience, it is suitable to use a spherical harmonic analysis, which 
is widely used when dealing with anisotropics on large angular scales: 

oo m = + / 

2£ 0/ (x 1 ) = E E A, m Y, m (n), (36) 

1=0 m = -l 

oo m = -\-l 
;=o m=-l 



and 

A lm = J 2£ /(xi)^* (fi^fi, (38) 

a lm = J 2E lf {£ 1 )Y* m {n)dn. (39) 

Substituting Eqs.(34) and (35) into Eqs.(38) and (39) respectively and using the Rayleigh 
equation 

/=oo 

e ,-k.x = e ^co S7 = i'(2l + l)j,(kx)P,(co S ^) (40) 



;=o 



and the addition expression 



4 m = + l 

Pi(cos^) = ^— £ Y lm (n k )Yf m (n), (41) 



where j\ and P\ are the l-th spherical Bessel function and the l-th Legendre function respec- 
tively, we obtain the following equations: 

Alm = ^I 8 ^ 2kH * r ) Y im(^)d^ (42) 



i 2l u 2 



air, 



f f |4 1 A k2 j,(2A; 1 ^ - 1 )j,(A; 2a; )^ ro (n kl )^ ro (n k2 ) 



12tt 5 

COsfcUX — id ■ x) ,o ,o, ,o, 

-J 3 xJ 3 k!J 3 k 2 . (43) 



According to random field theory, for a Gaussian random field, we have 

(S kl S k2 S k3 S u ) = | 8 kl | 2 | 4 3 | 2 <5(ki - k 2 )<5(k 3 - k 4 )+ | <5 fcl | 2 | 8 ki | 2 <5(ki - k 3 )S(k 2 - k 4 ) 

+ | 2 |4 2 r^kx-k^ka-k^. (44) 

Combining the above expansions and integrating for angular coordinates within the cone 
C(o?fi), we finally obtain 

<l a fm l> = <l a lmA^ m \> 


2i7 Q 2i7 Q oo oo 

+ 2 / / / / | 4 : | 2 | 4 2 | 2 jf(2k 1 HQ 1 )ji(k2X 1 )ji(k2X2)x 1 X2dx 1 dx2dk 1 dk 2 }, (45) 
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where the averaging is done over the entire sky and over all observation positions in the 
universe. Finally, the angular cross correlation function, generated by interaction between 
the gravitational lensing effect and the primordial perturbed part of the CBR, is as follows: 

oo 

°c(<x) = 7~ J2( 21 + !) 2 <l aim | 2 > Pi(a). (46) 
4?r 1=2 

For a double beam switch experiment, the observable 6T(a)/T } produced by gravitational 
lensing, is 

1 00 

6T{a)/T = — £(2Z + 1)2 <| a lm | 2 > 12 (1 - P,( a )). (47) 



IV. NUMERICAL RESULTS AND CONCLUSIONS 

We assume that the fluctuation spectrum is the Harrison-Zel'dovich spectrum when cal- 
culating the gravitational lensing effect on anisotropics of the CBR. For comparison with 
observation, an appropriate normalization is necessary. We use the rms mass fluctuation 
(SM/M) 2 = 1 within a sphere of radius r = 8h~ 1 Mpc as the normalization condition. Of 
course, how to select normalization conditions depends somewhat on how the mass distri- 
bution traces the galaxy distribution. At present, this is not clear and so the selection of 
the normalization condition may produce an uncertainty in the numerical results. However, 
this does not significantly affect our conclusion. With the above normalization condition, 
the normalized power spectrum is given by 

I 4 | 2 = r 4 k (48) 

In order to cancel the unknown nonlinear effect, we need to apply a low-pass filter function 
to truncate the spectrum. For mathematical convenience, we choose this as an exponential 
form e~ krt . Thus, for making calculations, the power spectrum is replaced by 

|4 \ 2 =r*ke- krt . (49) 

Here r t is a cut-off scale. Integrating Eq.(45) over k, we obtain 
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2H' 1 

( a lm) = o r I H- / Ql{ AT T-l x dx\ + 

o 7 

2_ff ( T 1 2_ff ( T 1 

^2 1 + r 2^2 o o r 2 + x 2 + ^2 

-z^— s - ^ y y g ' ( 4 2,,^) 2 ^ 2} 
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Wl 



ith 



/o = r /iLu = 3r " ( 51 ) 

J I 4 | d 3 k 

where Qi is the second kind of Legendre function. We have integrated Eq.(50) numerically 
in the case of r t = hh~ x Mpc } r = 9>h~ x Mpc and h = 0.75. The behavior of the multipoles as 
a function of the harmonic number 1 is showed in Fig.l. The predicted cross function of the 
temperature fluctuation of the CBR due to the gravitational lensing effect for comparison 
with double beam switch experiments is showed in Fig. 2. We have also made numerical 
calculations varying the parameters r t , r and h within acceptable regions. However, the 
numerical results are not sensitive to changes of these parameters. From Fig. 2, it is obvious 
that the effect of gravitatioal lensing on anisotropics of temperature of the CBR on the scale 
of arcminutes are too small to significantly amplify or depress the primordial anisotropics of 
the CBR. We can then safely conclude that the gravitational lensing effect on anisotropics 
of temperature of the CBR on scales of a few arcminutes are negligible. Our conclusion 
disagrees strongly with that of Cayon et al. who argue that gravitational lensing, produced 
by linear density perturbations, may enhances AT/T by ~ 20% or even more. 

In addition, the problem in the calculation presented by Cayon et al. (1993) could be 
that they considered the effect on the anisotropics of the CBR caused only by deflections 
of light lines, but not by convergences and divergences of light beams. In principle, the last 
effect could also produce a considerable change in the fluctuation of the intensity of the CBR 
compared with that caused by deflections of light lines. 
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Fig. 1: The hamonic amplitude of the secondary 
fluctuations produced by the gravitational lensing 
as a function of 1. x and H< represent the amplitude 
of the quadrupole moment of primary fluctuations 



for h=1.0 and h = 0.5 respectively (Holtzman 1989). 
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Fig. 2: The secondary temperature fluctuations as a function of 
angular, x and * represent the primary temperature fluctuations 
at 4.5' for h=1.0 and h = 0.5 respectively (Holtzman 1989). The 
arrows at 4.5', 7.2' and 30' show the upper limits of Uson & 



Wilkinson (1984), Readhead et al. (1989) and Gaier et al. (1992). 



